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Constant Proportion Portfolio Insurance (CPPI) is an investment strategy designed to give par- 
ticipation in the performance of a risky asset while protecting the invested capital. This protection 
is however not perfect and the gap risk must be quantified. CPPI strategics arc path-dcpcndcnt 
J>. and may have American exercise which makes their valuation complex. A naive description of the 

state of the portfolio would involve three or even four variables. In this paper we prove that the 
system can be described as a discrctc-timc Markov process in one single variable if the underlying 
^ asset follows a process with independent increments. This yields an efScient pricing scheme using 

transition probabilities. Our framework is flexible enough to handle most features of traded CPPIs 
including profit lock-in and other kinds of strategies with discrete-time reallocation. 
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1 Introduction 

Constant Proportion Portfolio Insurance is a dynamic strategy designed to give participation in risky 
assets while protecting the invested capital |Mer71] . This is achieved by periodically rebalancing 
between a risk- free asset (Zero-Coupon bond) and a risky asset (share, index, fund, fund of funds. . . ). 
In the simplest form, if the underlying asset has no jumps and if one can rebalance continuously, the 
final payoff depends in a deterministic way on the risky underlying. However both hypotheses are 
strong and do not fit real market conditions. If they are relaxed, the strategy is not as efficient: 
there is a small chance of not recovering at maturity the invested capital. This gap risk may be 
held by the issuer, so that the principal is really guaranteed to the investor. In this case, there 
is an option included in the product, which must be priced and hedged. For a very simple CPPI 
strategy with continuous rebalancing, the gap risk comes only from instantaneous jumps and can be 
quantified analytically |CT07j . With a discrete rebalancing scheme, there is a closed formula for the 
embedded option price if the underlying follows a Black-Scholes diffusion |BBM05j . This formula 
can as well be generalized to jump-diffusion models, and more generally to Levy processes |PL09j . 
However these methods work only for an idealized CPPI product where there are no caps on the 
risky exposure, no spreads on the risk-free and financing rates, no fees, no profit lock-in, a natural 
bond floor... Unfortunately, real CPPIs have usually such features, which prevent from using closed 
formulas. Moreover, such formulas does not hold for options with strikes differing from the guaranteed 
amount. This makes necessary to use other methods for real- life, more complex products. As CPPIs 
are very path dependent, they are usually priced through Monte-Carlo simulations (see [B K95| for an 
example) . Extreme value theory has also been used to estimate the gap risk of such products |BP02j . 

Monte-Carlo pricing is perfectly suited to path dependence but the dimension of the problem is 
generally quite large. As an example, a monthly CPPI defined on a single risky asset with monthly 
rebalancing and a 10 years maturity requires simulating 120 values of the underlying. This means a 
120-dimensional Monte-Carlo integration. Furthermore, the tails of the distributions are crucial in the 
pricing. First of all, the lower tail gives the gap risk part. It must be computed with enough precision 
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to produce a reasonable Put price and an accurate Delta. The upper tail is important to reach 
the correct mean value of the CPPI strategy. The payoff distribution is close to a shifted lognormal 
distribution with high leverage: most trajectories will end below the mean and will be compensated by 
a few very high terminal values. The high path dependence and barrier-like structure of the strategy 
makes the Put price and its Delta converge very slowly. As an illustration, we refer to figure [l] in 
which we present an example of the convergence of a put option on the CPPI portfolio. One needs 
100 millions Monte-Carlo simulations (and hours of computation) to get a rough estimation of the 
price of the derivative. 

An other way to price CPPI derivatives would be to use a PIDE scheme. This would also handle 
the American or more precisely Bermudan exercise which is more and more present in such products. 
However, there are at least three variables to propagate which depend on each other: the risky 
underlying value, the CPPI portfolio value and the risk asset weighting. If some profit lock-in feature 
is present, the guarantee level increases the dimension of the PIDE to four. The high dimensionality 
combined with the large number of time steps needed (at least the number of reallocations) makes 
this solution very computationally expensive. 

There are drastic simplifications if one considers the portfolio value at rebalancing times only. 
We prove that if the underlying asset has independent increments, the dynamics of the portfolio at 
rebalancing dates is described by a discrete time Markov process in one single variable. This property 
leads to a simple and efficient numerical scheme. The main advantages of this technique are its speed, 
the smoothness of greeks and the ability to handle Bermudan exercise. For example in the simplest 
cases, only a few milliseconds are needed to get an accurate price. More complex cases with profit 




10000 



100000 
Sampling Count 



le+06 



le+07 



le+Oi 



Figure 1: Convergence of the Monte Carlo method for the price of a put on a CPPI portfolio as 
functions of the number of paths. Maturity is 10 years, multiplier is 4: reallocation is performed 
monthly, volatility is 35% and initial investment is 1 000. The strike is set to this value. The dashed 



line indicates the best estimate obtained in section 3.4 
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lock-in, artificial cushion or coupons take more time but never more than a few seconds. 

In section [2j we start by describing our pricing method in simple cases for which the guarantee 
amount is constant: after introducing our working hypothesis, we formulate the main result concerning 
the transition probabilities of the CPPI portfolio at rebalancing dates. Then we discuss the numerical 
implementation of this method in section [3j In section |4j we allow the guaranteed amount to vary 
over time through a feature called "profit lock- in". Under some homogeneity assumption on the 
rebalancing rule, we show that the problem is still one-dimensional and discuss the steps involved in 
the numerical implementation. 



2 Basic algorithm 
2.1 Sketch of the method 

At a time t, we consider a portfolio (the CPPI portfolio) with total value C{t). This portfolio is split 
between a risky and a risk-free components: 

C{t) = C'^^^it) + C'"^'^-^'^<=^(t) 

The risk- free part is invested in zero-coupons of same maturity as the CPPI, B{t), whose dynamic is 
completely deterministic 

B{t') = 5(t)e''(*'-*) 

with r the risk- free interest rate (taken constant for simplicity). The risky part is invested in an asset 
S{t). For simplicity in this part we consider a Black-Scholes dynamics which under the risk-neutral 
measure is 

dS{t) = rS{t)dt + aS{t)dW{t) 

A CPPI strategy is designed to protect some part of the investment and we shall denote by G 
this "guaranteed" amount. At this point, G is just some constant, which is generally taken to be the 
initial value of the portfolio C(0). Taking a fraction 1/G of the portfolio, we can consider without 
loss of generality a unit guarantee. In this section we take therefore G = 1. 

Let to = < ti < ■ ■ ■ < tn he some fixed times, called rebalancing dates. At each rebalancing 
date ti, the allotment of the portfolio between the risky and risk-free components is adjusted (without 
any change in the total value of the portfolio) according to a rule specified through a non-negative 
deterministic function w{t,x) called the risky asset weighting (or simply RAW): 

G'''''y{ti) = w{u,c{u))c{ti) (1) 

For a CPPI strategy, w is usually chosen of the form 

(t,C7)=max(m^^^,o) (2) 



where m is called the multiplier of the strategy and H is a positive deterministic function called the 
floor whose terminal value is the guarantee (equal to 1 here). A natural choice is to set the floor to 
be 

H{t) = e-'^(*"-*) 

H{t) is equal to the value of a (unit nominal) zero-coupon bond one must hold to recover 1 at 
maturity tn- Within this setup, this investment strategy indeed offers some protection: when the 
portfolio poorly performs the investment is smoothly allocated to the risk-free asset. With continuous 
rebalancing and without jumps in the asset price, the bond floor cannot be breached. With discrete 
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rebalancing there is still a probability of going under the floor: the floor is breached if the risky asset 
goes down by more than 1/m between two rebalancing dates. 

Between rebalancing dates, the composition of the portfolio is left unchanged. For t G [ti,ti+i], 
this reads 



c{t) = c'''^y{ti 



'S'(^) , ^risk-free/. \ -^(^) 



(3) 



In the case t = tj+i, this equation gives the evolution of the portfolio between two consecutive 
rebalancing dates: 

C{ti+i 



S{U) 

or equivalently, eliminating C'^^^^y using equation ([T]), 



frisky ^.^.-^ S{ti+l) _|_ ^risk-free j-^. ^ 



B{ti 



C{U+i) = w{ti,C{ti))C{ti) 



+ {l-w{ti,C{U)))C{ti) 



B{t 



i+l) 



B{U) 



(4) 



Then, at date tj+i, the manager computes the new Risky Asset Weighting w{ti+i, C{ti^i)) and real- 
locates the resources. Once again, we emphasize that this rebalancing procedure has no instantaneous 
effect on the total value of the portfolio. 

In a first step, we look at the risk-neutral probability distribution of the portfolio value over a 
single rebalancing period. It can be seen from equation Q that C(ti+i) depends on C{ti) and on 
the ratio of the risky asset values 5(ti+i)/5(ti). However, if the process logS'(t) has independent 
increments (which clearly is the case of the lognormal dynamic we consider here), this ratio does not 
depend on S{ti) itself. If in addition we suppose that the risky asset has independent increments, the 
distribution of C(tj+i) conditionally on time ij depends only on C{ti): 



C iU+i) <y\C{ti) 

S{U+i) y-{l 
S{U) 



w{ti,C{U))) C{ti)B{ti+i)/B{ti) 



w{ti,C{U))C{ti 



C{ti) 



(5) 



This equation relates the distribution of C(ti+i) conditionally on C{ti) to the distribution of S (ti+i) / S (ti) . 
In the case of Black-Scholes diffusion that we consider, the cumulative function is 



S{t 



S{ti) 



< z 



'inz- (r- {ti+i 



(y^/ti 



+1 



(6) 



Combining these two equations we get an explicit formula for P C(tiJ^\) < y\J'ti as a function 

of C{ti). The process for the portfolio value at rebalancing times is therefore a Markov process in 
discrete time. The general framework for the study of Markov processes leads to introducing the 
risk-neutral probability density of C{tj) conditionally on C{ti), also called transition kernel for the 
period [ti,ti+i]: 

^ ^\c{ti+i) <y\C{ti) = 



dy 



which is computed explicitly using equations ^ and ([s]). 

Let us consider a European derivative product V written on the CPPI portfolio with terminal 
payoff at t„ given by 

V{tn,C{tn)) =P{C{tn)) 
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If we know the value of this derivative at time ti+i for ah values of C(ti+i), its value at time U and 
CPPI value C{ti) is given by the discounted risk- neutral expected value 



v{u,c{u)) 



g-r(ti+i-ti) ]g 



v{u+i, c{ti+i))\c{ti 



,-r{ti+i-t, 



Jdy ^taiU<^iti),y) V{ti+i,y) 



Starting from the terminal payoff and applying recursively this formula yields today fair price of the 
derivative V{to, C{to)). 

This backward recursion is also suited to Bermudan exercise (on rebalancing dates): 
F(ti,C(tO) =max(^e-'-(*»+i-*^) J dy ^taiUCiti),y) V{U+,,y) , P{C{U)) 

2.2 Mathematical formulation 

Wc collect the previous results (with slightly more general hypothesis) in the following proposition. 

Proposition 1. Let us consider a CPPI portfolio with a given set of rebalancing dates to,. . . ,tn and 
guaranteed amount G. If the following hypotheses hold: 

• the risky asset weighting is given by a deterministic function w{t,C/G); 

• the underlying risky asset follows a one- dimensional process, its logarithm has independent 
increments and it has a density 



(t>t,t'{z) = d. 

which depends only on t, t' and z; 

the risk free asset is deterministic with value B(t); 



S{t') 



[Sit] 



< z 



then the CPPI portfolio value taken at rebalancing dates C(ti) is a Markov process in discrete time. 
The probability transition operator of the normalized variable C{ti)/G 



-y 



G 



= X 



given by 



^UM+i{x,y) 



w 



{ti,x) 



X 



y - {I - w{ti,x)) xB{ti+i)/B{ti) 



w 



5{y -X 



B{t,+i) 
B{ti) 



if w{ti, x)x > 
if w(ti, x)x = 



(7) 



European derivatives on the CPPI portfolio are valued recursively through 

y(t,,C(t,))=e-'-(*^+-*^) /dy V{t,+,,Gy) 
and Bermudan derivatives using 

F(^^,C(^0) =max(^e-'^(**+i-*^) J dy ^t„u+,(^^,y^ V{ti+i,Gy) , P{C{ti)) 



(8) 



(9) 
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Proof. The proof goes along the hnes of the last section for a fraction 1 /G of the portfolio. The explicit 
expressions for the risk-free asset and the floor are simply relaxed and the lognormal cumulative law 
is replaced by a more general cumulative law. The Markovian property stems from the fact that the 
transition probability depends only on C{ti)/G. □ 

Remark 1. The condition on the density between rebalancing dates excludes stochastic volatility and 
more generally models with hidden variables. In fact such models can be used within our framework at 
the cost of increasing the dimensionality of the Markov process by the number of additional variables. 
For example a stochastic volatility model for the underlying will be converted into a 2-dimensional 
discrete-time Markov process. As matrix- vector multiplications on a grid of size have a complexity 
O(A^), adding a second dimension with grid size M will increase by a factor the computation 
time. Therefore it looks more reasonable to use a volatility regime-switching model rather than a 
model with continuous stochastic volatility. In such a model, there are a small number M of regimes, 
each with its own volatility. One computes the joint probability of the CPPI portfolio value and the 
volatility regime conditionally on the initial portfolio value and volatility regime. 



3 Numerical implementation 



3.1 Order 2 scheme 

The numerical implementation is rather direct from proposition [T] 

• The normalized value C{ti)/G of the portfolio was previously allowed to be any real number. 
This real line is discretized into a grid hj, j = 0...N. The (rescaled) initial value of the 
portfolio C{to)/G is supposed to lay on the grid node: C(to) = Ghj^. 

• The transition operator of equation ([T]) acting on real functions is replaced with a transition 
matrix M^^'*^"^^ acting on vectors. 

Let Qti,t^+i and Q[]^t^_^^ denote the cumulative density and partial expected value of the under- 
lying process 

Qu,u+Az) = IP 



' s{t') (s{t')y 
S(t) ^-°°''\s(t) J 



dy (f>u,u+iiy) 

dy (l)u,u+i{y)y 



Let Xj^ denote the value of the underlying return sit ) ^'-'^ which a portfolio worth Ghj at U 
is worth Ghk at tj+i: 

^(i) _ hk- {l- w{ti,hj)) hjB{ti+i)/B{ti) 
w[ti,hj)hj 

Starting from a grid point hj at time ti, for w{ti,hj^hj > the probability of being between 
hk and hk+i at time ij+i is the finite difference 



1) 



hk < 



G 



< h 



k+l 



G 
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Similarly the partial expected value is 



^i(v+i) 



E 

G 

w{ti, hj^hj 



G 



Cjti) 
G 



+ (l - w{ti, hj)) hj 



B{ti) 



For w(ti, hj^hj = we have respectively 



^jk 
^jk 



1 



Bjtj+i) 



B{k+i) 
' B{U) 



B{U) 



Using these matrices, we take as transition matrix between ti and tj+i 



M 



with 



M 



jk 



+ («,«+!) 
jk 



M. 



jk 



^- (M+1) 

jk jk 



h'k+iQjk - Hjk 
hk+i — hk 

^j,k-l 



hk — hk-1 



These matrices are such that the probability Qj*^*^^"* to be inside the interval [hk, hk+i[ and the 
partial mean Q 



i(j,j+i) 



jk 
hk+i 



are exact: 



hk+l 



dy ^(u,u+i)ih j,y) 



<iy^{u,u+i)(hj,y)y 



^jk 



For a European product, we define a vector Vj with the value of the derivative at date ti 
conditionally to C{ti) = Ghj. The backward propagation ([s]) is turned into the matrix-vector 
product 



M 



jk 



hk + Af 



- («,«+!) 



k 



with terminal condition 



The fair price of the derivative at time is given by 

Vito,G{to)) = vl^^ 
In the case of a Bermudan product the recursion is 

y«=max('e-'^(*»+-*»)5^M(;'^+^Vi*+^) , P{Gh,) 
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3.2 Error analysis 



We now sketch an analysis of the numerical error for a European payoff. Denoting by N the number 
of grid points and omitting the discount factor, the error in backward propagation on one time step 

from to ti is 



TV 

1 f\ *f 



k=0 

N-l 



M+(^'^+^) and M^: 

j^^^^ been defined such that they give exact results on constant and linear 
terms on each interval [hk] hk+i[. 



dy<^ti,ti+iihj,y) 



hk+1 — hk 



hk+i — hk 



As in the case of a trapezoidal quadrature, we are left with a second order term 

N-l 

12 



N-l p 



k=0 

with A/ifc = /ifc+i — hk- The second derivative can be evaluated by finite differences at the same order 
of precision. We denote it by D'^V^^\ 

The error computed at time ti and grid node hj must then be backward propagated to time to 
to get the contribution to the final price. Summing errors from all time steps and grid points we get 
the total error for a European payoff 



1 

12 



N-l 
fe=0 



.j=i 



As Ahk is proportional to AT ^, we end up with an asymptotic error 

eih,) = 0{N-') 

Moreover, minimizing the error under the constraint of a constant Ahk gives an optimal grid step 



Ahk 



1 



Y^e-'(''-"^MfJl^vt> 



i=l 



(10) 



i-1 



with M*^'^'*) the transition matrix between time to and tf. M^*^'*) = M'-'''^^'^^ In order to estimate 



1=0 



the optimal grid, one can either use an approximative analytical solution for the portfolio values and 
transition matrices or run a first round of the algorithm using a basic coarser grid as an initial guess. 
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3.3 Higher order schemes 



The algorithm can be generahzed to reach a 0{N~'p) convergence for any p. To construct the transition 
matrix in the second order case, we have weighted the boundaries of all grid intervals in order to 
match both the probability to be in the interval and the partial mean on the interval. Introducing 
p — 2 intermediate grid points inside all intervals allows to match all conditional moments on every 
interval up to order p — 1. Following the lines of the error analysis of section 3.2 one gets an 0{N^p) 
error 



e oc 



N-l 

E 

fc=0 



with an optimal grid for a European payoff given by 



A/ifc oc 



[ID 



There is however a drawback: this (partial) moment matching procedure may introduce negative 
transition probabilities on some grid points, which could result in numerical instabilities. 



3.4 Numerical results 

We consider a put on a CPPI portfolio with strike at the guaranteed amount. Maturity is ten years, 
the portfolio is reallocated monthly with multiplier 4; the bond floor is linear starting at 75% so that 
the initial investment in the risky asset is 100%. The initial portfolio value is 1000. On rebalancing 
dates, proportional fees are removed from the portfolio value at an annual rate of 30 bp. For our 
test we choose a constant, flat 35% volatility. The yield curve is the EUR yield curve on January 15, 
2010. 

Prices and computation times (CPU time) as functions of the grid size are presented in table [l] 
and figure [2] for four numerical schemes: second (II) and third (III) order schemes using the grid 
described in the next paragraph and second (IF) and third (IIF) order schemes using the optimal 



grids computed using formulas (10) and (11) with p = 3, bootstrapped from the previous grid. 

The grid used in cases (II) and (III) is built in the following way. We first evaluate an upper 
bound inverting the cumulative distribution of a self-financing CPPI with natural bond floor and 
continuous reallocation: this distribution is a shifted lognormal distribution with volatility ma. We 
set the bound such that the cumulative at this bound is 1 — 10"^''. In the worst case where the risky 
asset drops to zero, the loss is equal to m times the current portfolio value. Hence we set the lower 
bound to (1 — m) times the upper bound. These bounds are usually very large but the grid must 
be refined where there is convexity in the CPPI price, according to the error analysis of the previous 
section. We therefore take a parametric form for the grid between the bounds. We split the grid in 
three parts: a log-linear (y = ae^^) central region which includes the bond floor value at all times and 
the strike of the option or the guarantee and two log-quadratic [y = c + de^^ ) outer regions. Regions 
are glued together so that the parametric function describing the grid is C^. This particular choice 
of the parametric form allows us to have sufficient refinement in the central region while reaching 
extreme bounds. 

We observe that for all schemes, prices smoothly converge as the grid size is increased. As expected, 
third order schemes perform better than second order ones. Asymptotically, the computation of an 
optimal grid gives better results. However at lower grid sizes the error is larger; in addition the 
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Figure 2: Price of a put on a CPPI portfolio as a function of the grid size in four implementations of 
the algorithm: (II) second order scheme, (W) second order scheme with computation of an optimal 
grid, (III) third order scheme, (IIP) third order scheme with computation of an optimal grid. 
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Table 1: Price of a put on a CPPI portfolio (top numbers) and CPU time in seconds (bottom) 
as a function of the grid size in four implementations of the algorithm: (II) second order scheme, 
(IF) second order scheme with computation of an optimal grid, (III) third order scheme, (IIP) third 
order scheme with computation of an optimal grid. 

computation is almost three times longer. In practice, we therefore prefer scheme (III) which gives a 
smooth and fast convergence. 

A precise price can be obtainec|^in less than Is. When the rates or the volatility are not constant, 
one needs to compute a transition matrix for each rebalancing period. In our test case, this means 
we generate 120 matrices. When all parameters are constant or piecewise constant, the number of 
matrices to be computed can be greatly reduced and the pricing is faster. 

4 Profit lock-in 

4.1 Pricing with profit lock-in 

CPPI strategies often include a profit lock-in feature: a fraction A of the performance is periodically 
locked in. G therefore becomes a piecewise constant stochastic process G{t). The floor is scaled 
accordingly as 

H{ti) = e-<'--'^^G{ti) 

such that the terminal value of the CPPI portfolio should ideally be at least G{tn)- 

We denote by , / = 1 . . . the subset of rebalancing dates at which the profit lock-in occurs 
(typically annually). For convenience, we set ^Q = Q and = n. Whenever the performance of the 
portfolio is positive between t^j and t^j^-^ the guarantee is raised to 

G{te,^,) = G{te,) + A [C(t,,^J - C(tf,)]+ (12) 

^Computations were performed on a standard PC equipped with an Intel Core 2 6600 at 2.4 GHz and 3.50 GB of 
memory. We used the optimized Intel's Math Kernel Library to perform matrix-vector multiplications. Quoted times 
are CPU times. 
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This expression can be recast into the more general form 



with f{x,y) a positive function. 
The rebalancing rule 



is of the more general form 

^""'('■) 

We are concerned with the pricing of a European derivative product V written on the CPPI 
portfolio whose terminal payoff P might depend on C and G at maturity. The fair price at time ti^ 
can be written as the risk-neutral expectation 

V{te„C{te,),G{te,)) = e~'-(*^^+i-*^^) E [v{te,^„C{te,^,), | Tt,^] (14) 

This problem is two-dimensional as one must know the joint law of C and G to apply this formula. 
However, beautiful simplifications occur in the case of a homogeneous payoff of the form 

P(C,G) = GP(^) (15) 

This includes the CPPI portfolio itself or options with strike in percentage of the final guarantee. 
Owing to the special form of the payoff we can state the following proposition. 

Proposition 2. Consider the same hypothesis as in proposition^ except for the fact that the guaran- 
tee is not constant but is rescaled at dates t£j along equation ( [13^ . For a European derivative product 
whose payoff at maturity is given by equation (15), the fair price at date tij is homogeneous: 

y(t,„C(t,J,G(t,,)) =G(t,,)y(t,„g|^) (16) 
and V satisfies the recursion formula 

with the kernel 

(x, y) = j Az5{y- fi^, ^) z) (18) 

In this equation, ^tej,tij^^{x, z) is the transition operator defined in propositionui 



Proof. At time i„, equation (16) is verified due to equation (15). We assume now equation (16) 
holds at time tij_^-^^. We denote by tj^^^ a time just before the profit lock-in, i.e. before the jump in 
G{t), but after possible jumps in the risky asset S{t) (this describes the state of the system before 
applying the lock-in). As there is no profit lock-in between t£j and i^^^, we can use proposition jlj 
with transition operator ,- and write 

V{te„C{te,),G{te,))=e-^^'^''^^-'^^^ J dy v{ti^^^,G{te,)y,G{te,)) (19) 
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When the lock-in is performed, the portfoho value is not changed: 



This value can be rewritten further using the recursion hypothesis (16) as 



V(tJ^^^,G{te,)y,G{te,)) = G(vJ F V, 



(20) 



The profit lock-in is deterministic as a function of G{t£j), C{tij) and C{t£j_^-^). It is given by 



equation ( 13 ) which reads here 



G( 



(21) 



Using equations (20) and (21), equation (19) is rewritten as 



V{te„Gite^),G{te^)) = 



y 



Inserting 



dz 5 \ z 



y 



f(^ V 



we find that C(t£^), is indeed of the form given by equation (16) with equations (17) 

and (18) satisfied. □ 



Remark 2. Proposition [2] can be generalized to a homogeneous payoff of degree a: 

P{G,G) = G''p{^ 



The only modifications are in equation (18) where G{tij) is replaced with G{tej)°' and in equation (18) 
where f{x, z) is replaced with f{x, z)". 

Remark 3. Proposition [2] can also be generalized to a Bermudan payoff with exercise allowed at lock- 
in or reallocation dates: as usual, the maximum of values with and without exercising has to be taken 
during the backward propagation. 

4.2 Numerical implementation 

The numerical implementation of proposition [2] is straightforward, building on the basis of section [3| 

• A grid is built as in the basic case of section [sj such that C(to)/G'(to) = hjg for some index jo- 

• Prom the local transition matrices Mjj^ of section 1, we compute the transition matrices 
between two lock-in dates as the matrix product 
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Using the lock-in function f{x, y), the discrete analogue of the kernel ^tij,t(j^^ (x, y) of equation 



(18) is a matrix Mjj^'^~^^\ We compute first 



J hk 



dzh 



hk+i 



''k 



We define the matrix M^^'-^^^^ from Q^^'^^^^ and Q^^^'^~^^^ exactly as in section p 



with 



~+ (1,1+1) ^ n^k+i^jk ^jk 

hk+i - hk 

Mjk = — 1 1 



IS discretized as vj^^ = V{tij,hj). It is propagated 



• The rescaled payoff V of proposition 
backward from lock-in date to lock-in 3ate through equation (17) which becomes 

yd) ^ y-M(V+i)^(/+i) 

J / ^ jk k 

k 

with terminal condition 

• The fair price of the derivative at time is given by 

V{to,Cito),Gito)) = Gito)vi;^^ 



4.3 Probability transition 

Proposition [2] provides a pricing algorithm for payoff which are homogeneous in C, G. The distribution 
of the portfolio value can also be computed, which may be useful for valuation and risk management 
of options with strike given as a fixed amount or for computing the final probability density function 
of a CPPI portfolio. 
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Proposition 3. The probability transition operator between two lock-in dates tgj and 



C{t, 



G{ti>j.) 



can be computed from the transition operator without profit lock-in as 



t,,tiAx,y)= [ dz ^ $t t (x,2:) ( A (22) 

' ^ J f[x,z) 'I'^i+i 'w'J\f{x,z) f{x,z)J 



Proof. This result is a rewriting of 



E 



y 



dz E 



using equation (flSl). 



X = 

CM. 

G{h^) 



E 



^, C{t,,) Gjt,^^,) ^ 
G{tej_^_^) G{tij) 



G{tei+i) G{tij^^) 
G{tej_^_^) G{tij) 



□ 



Supposing for simplicity that 
tn can then be obtained as 



E 



6{G{tn)-y) I C{to 



and £]sr = n, the distribution of the portfolio at final maturity 
1 



G{to) 



to, in 



X y 
GitoYGito] 



(23) 



The main difference with the homogeneous payoff is the following: whereas the dynamics of 
G{t£j)/G{tij) taken at lock-in dates is Markovian, it is not true for the dynamics of C{t£j) itself. 

The numerical implementation is similar to the other cases. Equation (22) is not a product of 
operators. As a consequence, it is also numerically more time-consuming. 



5 Numerical examples 

In order to illustrate the smoothness of the results, we plot in figure [3] the probability distribution at 
maturity for three different CPPI strategies. We take a risky asset which follows the jump- diffusion 
model introduced by Kou |Kou02j . There are upward and downward jumps, both exponentially 
distributed with distinct parameters. The diffusion has volatility 20%, downward jump with intensity 
0.1 yr~^ and mean value 0.1, upward jumps with intensity 0.1 yr~^ and mean value 0.05. We consider 
strategies on this risky asset with a maturity of 10 years and weekly rebalancing. The risk-free rate 
gives an initial zero-coupon value of 0.606. 

• We consider first the vanilla case of section [2] with multiplier m = 4: 

/ f C-H(t) \ 
w[t, 6 j = maxi m — , 1 

• We add a cap (150%) to the risky asset exposure so that 

f ^^ ■ ( r. ( G-H(t) \\ 
[t, G j = mm( Cap , maxI m — ,0 1 I 



• Removing the cap, we finally incorporate to the strategy an annual lock-in of 75% as described 
in section m 

We refer the reader to |PL09| for further discussions of the effects of the various features commonly 
found in CPPI-based contracts. 
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Figure 3: Probability Density of a CPPI portfolio value at maturity for three different strategies. 
For the first curve, the CPPI has maturity 10 years, weekly rebalancing, multiplier 4- The initial 
investment and the initial guaranteed level are 1. The mean price of the strategy is 1.65. The risky 
asset follows a Brownian motion with 20% constant volatility with additional jumps (Kou model): 
downward jump with intensity 0.1 yr""*^ and mean value 0.1, upward jumps with intensity 0.1 yr~^ 
and mean value 0.05. For the second curve the risky asset weighting is capped to 150%. For the 
third one there is no cap on the exposure but a profit lock-in of 75% of the performance is performed 
annually. The second graph is a zoom on the tails of the distribution. 
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6 Conclusion 



We introduced an efficient scfieme for tfie pricing of CPPI strategies and options. Instead of following 
the underlying spot, the portfolio value, the risky asset weighting and the guaranteed amount as in 
a Monte-Carlo pricing or a classical PIDE backward propagation, only the CPPI portfolio value at 
rebalancing dates is considered (or the ratio of the CPPI over the guaranteed amount in the case 
of lock- in). We proved that under the hypothesis of independent increments of the logarithm of 
the underlying, the portfolio value at reallocation dates (or lock-in dates) is a discrete time Markov 
process in one single variable. We derived a natural pricing scheme which uses this property to price 
both European and Bermudan derivatives on the CPPI portfolio. Numerical experiments exhibit a 
very fast convergence. This can be used to estimate the gap risk of guaranteed CPPI products, price 
options on CPPI strategies and design hedging strategies. 
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